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Detection of two-photon exclusive production of supersymmetric pairs at 
the LHC 

Nicolas SchuP * and Krzysztof Piotrzkowski a 

a Universite catholique de Louvain, Center for Particle Physics and Phenomenology (CP3), 
Louvain-la-Neuve, Belgium 

The detection of pairs of sleptons, charginos and charged higgs bosons produced via photon-photon fusion at the 
LHC is studied, assuming a couple of benchmark points of the MSSM model. Due to low cross sections, it requires 
large integrated luminosity, but thanks to the striking signature of these exclusive processes the backgrounds are 
low, and are well known. Very forward proton detectors can be used to measure the photon energies, allowing 
for direct determination of masses of the lightest SUSY particle, of selectrons and smuons with a few GeV 
resolution. Finally, the detection and mass measurement of quasi-stable particles predicted by the so-called sweet 
spot supersymmetry is discussed. 
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1. The two-photon exclusive pair produc- 
tion and search for new physics 

The 77 production of pairs of charged massive 
particles offers interesting potential for searches 
of the beyond standard model (BSM) particles. 
In a recent paper pQ, the initial comprehensive 
studies of high energy photon interactions at the 
LHC were reported. In the present contribution, 
the selected results discussed in ref. Q] are in- 
troduced and supplemented by new results. The 
exclusive two-photon production, pp — > pXp, 
provides clean experimental conditions, thanks 
to absence of the proton remnants. Well defined 
final states can be then selected, and precisely 
reconstructed. Moreover, detection of the two fi- 
nal state protons, scattered at almost zero-degree 
angle, in the dedicated very forward detectors 
(VFDs), provides another striking signature, ef- 
fective also at high luminosity and with large 
event pile-up 2 3\. In addition, the photon ener- 
gies can be then measured and used for the event 
kinematics recontruction (see section I3.2|l . 

The two-photon cross sections of pair pro- 
duction in general are determined by the mass, 
spin and charge of the produced particles, so 



the rate of produced particles at the LHC can 
be well predicted using the Equivalent Photon 
Approximation for the photon fluxes [4]. One 
needs however to apply corrections due to a pos- 
sibility of strong interactions between protons, 
the so-called rescattering effects. The resulting 
suppression of the cross sections depends in par- 
ticular on the invariant mass of the exclusively 
produced state X, and for the processes discussed 
in the following, as pp —> pWWp, it is estimated 
to be about 10% [5]. This correction is ignored 
in the present analysis. In table [TJ the two- 
photon production cross sections computed using 
MadGraph/MadEvent [6] for pairs of W bosons 
(spin 1), and of fermions (spin 1/2) and scalars 
(spin 0) for a couple of masses are shown [T]. 
One should note that if the inelastic production 
is considered, when one of the protons dissociates 
into a low mass state, the cross-sections increase 
by about a factor of three [2] . 

Supersymmetry predicts new scalar and 
fcrmionic particles above the W mass scale, 
and is therefore a prefect candidate for such 
novel and complementary searches using photon- 
photon collisions at the LHC. Two photon pro- 
duction of non strongly-interacting SUSY par- 
ticles at the LHC was discussed in 7 . In the 
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Table 1 

Cross sections for several examples of the exclu- 
sive two-photon pair production at the LHC. (F for 
fermion, S for scalar). [I] 



Produced pairs 


mass [GeV] 


a [fb] 


w+w- 


80 


108.5 


F+F- 


100 


4.064 


F+F- 


200 


0.399 


s+s- 


100 


0.680 


s+s- 


200 


0.069 



present work, the experimental aspects of detec- 
tion of such events at the LHC are studied. One 
should stress here, that in general the diffractive 
exclusive production of non strongly-interacting 
particles is heavily suppressed and can be safely 
nee lecteoll. 

2. Detection of supersymmetric pairs 

In two-photon processes, in general the produc- 
tion and decays mechanisms are simple, there- 
fore not involving cascade decay problems present 
in many SUSY studies at the LHC. In order 
to ensure the cleanest event signature only fully 
leptonic final states are analyzed in the follow- 
ing. In the Minimal Supersymmetric Standard 
Model (MSSM) the lightest supersymmetric par- 
ticle (LSP) is the stable neutralino Xi- As it is 
shown in figure[l] the final state can be then com- 
posed of a pair of leptons and at least two neu- 
tralinos which escape the detection, together with 
two neutrinos for some decays modes. 

The requested final state for this MSSM search 
involves then: 

• 2 leptons (ee or e/z or /z/n) of opposite charge 
detected in the central region, 

• 2 scattered protons detected in the VFDs 
(assuming the tagged photon energy range: 
20-900 GeV), 



2 For example, the gluon mediated exclusive production of 
W boson pairs, is about 100 times smaller than the two- 
photon production at the LHC [5]. 
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Figure 1. Most relevant Feynman diagrams for 
the MSSM analysis. Left: Chargino decay into 
a charged (neutral) scalar and a neutral(charged) 
fermion. Right: Slepton decay into a charged lepton 
and a LSP. 



• missing energy due to escape of neutrinos 
and neutralinos, 

• acoplanarity. 

The only irreducible background process for this 
event topology is the exclusive two-photon pro- 
duction of pairs of W bosons, followed by leptonic 
decays. Indeed, e + e~, /-t + /i~ and t + t~ pairs pro- 
duced via 77 fusion can be easily supressed by 
requiring large acoplanarity and/or large missing 
energy. The cross section for the two-photon pro- 
duced WW pairs decaying leptonically reaches 
almost 7 fb. 

In sections [3] and [4j the analysis is presented 
of two low mass benchmark points in the usual 
MSSM plane (LM1 and LM9) [8], chosen to study 
the potential for detection of SUSY particles and 
of determination of their masses. For LM1 and 
LM9, respectively, the slepton and chargino con- 
tributions are dominant. Third scenario dis- 
cussed in section [5] involves the f-j 4 as the NLSP, 
and the analysis method described above does not 
apply. Signal as well as background samples were 
generated using the modified CalcHEP [9j and 
then passed to the modified Pythia [TO] , where the 
decay and hadronisation steps were performed. In 
order to simulate acceptance regions of the LHC 
central detectors, cuts (UJ) were applied on leptons 
at generator level: 

p T (e ± ) > 10 GeV, p T (^ ± ) > 7 GeV 
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1*7(^)1 < 2-5. (1) 

3. LM1 benchmark scenario 

Sparticle masses in the LM1 scenario are given 
in tabled] Application of the acceptance cuts (1) 
reduces the cross sections by factors between 1.5 
and 50, depending on the particle type. 



Table 2 

Masses of LM1 MSSM particles derived from running 
the RGE for m = 60 GeV, m 1/2 = 250 GeV, tg(/3) 
= 10, A = 0, fi > 0. (£ = e,n) 



mass [GeV] 





118 


xt 


178 


It 


187 


xt 


360 


ft 


111 




381 


ft 


190 


X? 


96 



Because the branching ratio BR(£ + — > £ + + Xi) — 
100%, selectron and smuon pairs are the ma- 
jor expected contribution in the analysis of di- 
leptonic final states. In addition, almost 65% of 
right sleptons and 75% of left ones fall within the 
acceptance window (JTJ. On the contrary, pairs 
of charginos (left diagram of figure [TJ) have other 
significant decay modes (including hadronic de- 
cays) therefore only 7% of the produced chargino 
pairs will be detected. Staus, since they decay 
into tau leptons, will produce mostly two r-jet 
final states, hence only 1 and 6% of the fi and the 
T2 pairs, pass acceptance cuts, respectively. One 
could also consider final states with one lepton 
and one r-jet or with two r-jets, but the statis- 
tical improvement is then only relevant for large 
tg(/3) models, where couplings to (s)taus are en- 
hanced. Finally, no relevant signal of the charged 
higgs boson pairs can be seen, since H + — > bt is 
the dominant decay mode. 

3.1. Lepton flavour sharing 

By considering only the information from the 
central detector and applying acceptance cuts, 



one already reaches a signal-to background (S/B) 
ratio close to 1/5, see table [3] This ratio can be 
easily improved by considering only same flavour 
dileptonic events. Indeed, because the signal is 
dominated by in and £l pair decays, around 
90% of the LM1 events is composed of e + e~ and 
leptons (see right diagram of figure [lj. At 
the same time the background W pairs are sup- 
pressed by a factor 2 by selecting same flavour 
leptons. 

3.2. Very forward detector information 

In contrast to the nominal pp studies which 
can only use kinematical quantities measured 
with the central tracker, calorimeters, etc., the 
main experimental advantage of the two-photon 
processes relies on the detection in VFDs of the 
two forward scattered protons [2|llj . To eval- 
uate impact of these detectors on the analysis, 
a complete set of VFDs, situated at 220 m and 
420 m on each side of the interaction point, is 
assumed. In this case, the accepted or tagged 
photon energy is contained between 20 and 900 
GeV. The efficiency of detecting both forward 
protons is then around 80% for the LM1 bench- 
mark events, so almost whole relevant photon 
spectrum is probed. Finally, one should note 
that the proposed forward detectors are capable 
of running at high luminosity allowing for se- 
lection of such exclusive dilepton events also in 
presence of large event pileup [3] . 

In order to simulate the photon energy resolu- 
tion, as can be obtained using the VFDs [12], a 
gaussian smearing of the photon energy with the 
following simple parametrization of its standard 
deviation: 

cr^ = max(1.5 GeV, £ 7 /100). (2) 

Useful quantities used to discriminate signal and 
background are listed below. An obvious vari- 
able which can be computed is the total photon- 
photon center-of-mass energy, 

W 11 = 2^jE 11 E 12 (3) 

where E lx , E 12 are the reconstructed energies of 
two colliding photons. The expected W 77 dis- 
tribution is shown in figure [2] for an integrated 
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where En, Ei 2 are the measured leptons energies. 
A conservative correction is made to account for 
the bremsstrahlung in electronic decays. It is 
assumed that the soft bremssthralung photons, 
Pt{i) < 10 GeV, are not detected. This results 
in the biased E m i SS in a small fraction of events, 
but otherwise leptons are very well reconstructed 
in the central detectors. Therefore, it is assumed 
that the energy and momentum of the leptons 
are known exactly, and the resolution of the re- 
constructed kincmatical variables is dominated 
by the photon energy resolutions. 



Figure 2. Distribution of two-photon invariant mass 
W 11 for the LM1 benchmark and integrated lumi- 
nosity L = 100 fb _1 . Two visible peaks are due to 
production thresholds of ij^jt^ and pairs. Veri- 

ous contribution are added cumulatively. The back- 
ground distribution of WW pairs is shown separately, 
and is rescaled to obtain similar size as signal. 



luminosity of 100 fb _1 . One can see two signif- 
icant peaks due to the production thresholds of 
right slepton (around 250 GeV) and of left slep- 
tons (around 400 GeV). In this way, the slepton 
mass spectrum can be probed by measuring the 
threshold energy for each peak, which is approx- 
imately equal to the sum of masses of the two 
produced sparticles. It should be stressed that 
the mass determination in this method depends 
only on the VFD energy resolution, and not on 
the resolutions of the central detectors. Actual 
precision of and 1^ mass determination is then 
mostly driven by available statistics. Moreover, 
this quantity can also be used to suppress the 
background since the IV 77 distribution for W 
pairs is well known and starts at about 2mw- 
The choice of the analysis cut W™ m could be 
then changed according to the SUSY mass con- 
straints coming from pp studies. 

The missing energy carried away by the neutri- 
nos and the ncutralinos can be estimated as 



The missing invariant mass distribution can 
then be defined as 



W„ 



= \ E. 



P 2 

miss 



(5) 



where P m i SS is the event missing momentum and 
is calculated in analogy to E m i SS . Missing mass is 
on average larger for the SUSY event sample since 
a supersymmetric event will always produce at 
least two massive LSPs. This can be seen in figure 
[3] where the expected distributions of W m i SS for 
100 fb _1 is shown. SUSY distribution is peaked 
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Figure 3. Distribution of missing invariant mass 
W m iss for the LM1 MSSM benchmark for the inte- 
grated luminosity L = 100 fb _1 . It starts at about 
2 itilsp for SUSY, at zero for the WW background. 



Pjniss — -^71 ~h -^72 



(4) 



slightly above 200 GeV, which is twice the mass 
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of the lightest neutralino in this model, while it 
starts at zero for the SM contribution (which is 
not rescaled here). This quantity allows for mea- 
suring the LSP mass with a high resolution and 
to suppress the W + W~ background by requiring 
a cut. 

3.3. Estimate of the LM1 signal signifi- 
cance 

To reconstruct W 77 and W m i SS quantities, it is 
assumed that the both forward scattered protons 
are detected in the VFDs within the full tagging 
range. Then, by applying the following analysis 
cuts 

• acceptance cuts (Q}, 

• W-y-y > 200 GeV and W miss > 186 GeV, 

• lepton angular cuts, Arj < 2.1 and AR < 
3.0 (R is the distance between two leptons 
on the pseudorapidity-phi angle plane), 

• same flavour lepton selection, 

one reaches a S/B ratio close to 2. The various 
contributions for the signal are given in table [3] 
and shown in figure |4] for the integrated luminos- 
ity L = 100 fb"\ and S ~ 51 and B ~ 26 events. 



Table 3 

LM1 signal and WW background cross sections be- 
fore and after the acceptance cuts (including the fla- 
vor selection), and after the analysis cuts. Values are 
given in fb. (I = e, fi. i= 1,2). 



Processes 





&acc 


&acc-\-ana 


77 "> t + J R 


0.798 


0.522 


0.403 


11 -> itfi 


0.183 


0.135 


0.089 


11 rffr 


0.604 


0.054 


0.003 


11 -> XiXi 


0.642 


0.043 


0.014 


77 -> H+H- 


0.004 


/ 


/ 


77 -► w+w- 


108.5 


3.820 


0.255 
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Figure 4. Missing energy distribution for the LM1 
benchmark and the WW background for the inte- 
grated luminosity L = 100 fb - . 



The 5<T discovery for the LM1 left and right 
sleptons is then reached already after 25 fb -1 
thanks to strong suppression of the irreducible 
background. It could still be improved by using 
additional cuts exploring the correlation between 
W 77 and Wmiss as it is done for the LM9 study 
below. Finally it could be improved even further 
by including the inelastic two-photon production, 
in this case however only one proton is detected 
and the kinematical reconstruction is not so effec- 
tive. For the same benchmark point, the nominal 
proton-proton studies claim 5a discovery after 
about 10 fb -1 [13j. However, determination of 
sparticle masses in this case is much more com- 
plicated. 

Similar two-photon analyses can be done for 
other benchmark points with low slepton masses 
as LM2, LM4 and LM6. 

3.4. Mass measurement 

The main advantage of the two-photon analy- 
sis is large sensitivity to sparticle masses. Mass 
determination using the production threshold 
values in W / 77 and W m i SS distributions is lim- 
ited by the number of selected events. Another 
approach, based on other kinematic quantities, 
can provide a method to measure mass of the 
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sleptons more on the event-by-event basis. 

The two-dimensional plots at figure [5] represent 
event distributions on the W 77 -, W m i SS plane for 
the MSSM processes and for the W + W~ back- 
ground, after acceptance cuts. One can observe 
that for the signal events these two variables are 
strongly correlated, and much less so for the back- 
ground. Moreover, it can be shown that the width 
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Figure 5. Scatter plots for the LM1 signal and the 
WW background on the Wmiss, W 77 plane. Up- 
per: AflAij events (orange/light) and W + W~ events 
(red/dark). Lower: p^fi£ events (green/light) and 
W + W~ events (red/dark). 



of the distribution is related to the mass of the 
produced sparticles. This demonstrates a close 
relationship between the MSSM masses, the 77 
invariant mass and the missing mass. An empir- 
ical quantity has been built in order to take into 
account this observation: 

{2m reco f = - {[Wl iss - 4m| ; ] V2 

+[^-H/ /2 ) 2 

where Wi ep is the invariant mass of the two lep- 
ton system, and mi ep is the lepton mass, and 
2m reco is the reconstruted mass of the produced 
sparticles. That relation does not work well for 
f and x pairs because in general they decay into 
final states with more neutrinos and neutralinos. 
The only needed input in this method is the value 
of the LSP mass, which can be derived from the 
Wmiss distribution. 

The reconstruction power of this empirical 
quantity is illustrated in figure [6] for the inte- 
grated luminosity L = 100 fb _1 . A narrow peak 
centered at 2m reco = 236 GeV — 2 x 118 GeV, al- 
lows for efficient and direct determination of the 
and jx^ mass. A second peak, centered at 
2m reco = 370 GeV = 2 x 187 GeV, with larger 
width, corresponds to and p, L pairs but is not 
so well visible. Right selectron and smuon mass 
might be determined using this empirical method 
with a few GeV resolution. In summary, this 
techniques of reconstruction of sparticles masses 
are very promising for the low SUSY mass scenar- 
ios, and it should be further studied using more 
complete, or even full simulations of the central 
detector performance. 

4. LM9 benchmark scenario 

In this second scenario, photon-photon pro- 
cesses would probe the MSSM signal via the de- 
tection of light chargino pairs. In LM9, in order to 
keep the mass for the lightest neutral higgs above 
the LEP limit, scalar sparticles are very massive 
(around 1.5 TeV, see table [4j> and consequently 
can not be detected using the two-photon events. 
The only significant contribution for the LM9 
signal arises from the lightest chargino detection. 
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Figure 6. Cumulative distributions of the recon- 
structed mass 2m reco for the LM1 signal and the WW 
background for the intergrated luminosity L = 100 
fb" 1 . 



W 77 > 220 GeV and W miss > 180 GeV, 



• W 77 < 1.5 W„ 



• Wi ep < 170 GeV, 

• no flavour selection, 

• angular cuts AR <2.6 and An <3.5 (R 
is the distance between two leptons in the 
pseudorapidity-phi angle plane), 

Additional cuts, in comparison to the LM1 study, 
the cut on two-lepton invariant mass and the cut 
using the correlation between W 77 and W m i SS 
are applied in order to suppress two-photon WW 
events, in situation where the flavor selection is 
not effective. The expected visible cross sections 



Table 4 

Masses of the LM9 MSSM particles derived from run- 
ning the RGE for m = 1450 GeV, m 1/2 = 175 GeV, 
tg(/3) = 50, A = 0, > 0. (£= e,n) 



mass 


GeV] 


it 


1450 


xt 


107 


*t 


1450 


xt 


223 


ft 


1054 




495 


ft 


1267 


X? 


65 



> 
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In that particular scenario, the branching ratio 
BRixt -> Xi + W*+) = 100%, and the signal 
and background decay schema are almost the 
same, with additional missing energy in SUSY 
events. After applying only the acceptance cuts 
(JTJ on generator level, one finds for the LM9 point 
a S/B ratio of about 2.5% with about 10 xtxi 
events for the integrated luminosity L = 100 fb -1 
(see first columns of table [5J . 

As for the LM1 study, one can also apply anal- 
ysis cuts on W 77 and W m i SS quantities as it can 
be deduced from figure [71 Best values were found 



Figure 7. Scatter plots for the LM9 signal and WW 
background events on the Wmiss, W 77 plane. 



after applying these analysis cuts are quoted in 
table [3 and result in respectively 7 and 41 events 
for the MSSM signal and the WW background for 
100 fb -1 . Those stringent cuts improve the S/B 
ratio to above 15%, without decreasing dramati- 
cally the number of signal events. Another analy- 
sis method taking into account semi-leptonic and 
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Table 5 

Cross sections before and after the acceptance cuts, 
and after the analysis cuts for the LM9 signal and the 
WW background. Values are given in fb. 



Processes 


o 


&acc 


&acc-\-ana 


77 -> X^Xi 


3.281 


0.099 


0.071 


77 -> xtx2 


0.273 


0.001 


/ 


77 W+W- 


108.5 


3.804 


0.408 



protons. For example, the two-photon invariant 
mass W 71 as shown in figure [8] for the integrated 
luminosity L = 100 fb -1 , assuming detection if 
both staus with pr{j\ ) > 10 GeV, and of two 
forward protons in the full VFD setup. The stau 
mass can then be directly measured by compar- 
ing the two-photon invariant mass with the stau 
momenta. Such a event-by-event mass measure- 
ment has a very good resolution, of about 5 GeV. 
Finally, the stau spin could be determined by 



fully-hadronic final states should then be consid- 
ered in order to tackle these particular supersym- 
mctric scenarios involving chargino pairs as dom- 
inant contribution. Consequently, the event sig- 
nature involving jets and r-jets would be more 
difficult to extract. 

5. Sweet Spot SUSY scenario 

In this last scenario, SUSY search strategy 
would be changed due to different phenomenol- 
ogy. In some non-common SUSY theories, the 
LSP is not Xi as m MSSM but rather the grav- 
itino. One interesting framework in which that 
occurs is the Sweet Spot SUSY [14], where the 
next-to-LSP is predicted to be the lightest tau 

, with low mass around 116 GeV. It would 
be quasi-stable, with a decay time of O(3000s). 
Production of pairs in photon-photon inter- 
actions will then be detected by observing very 
exclusive final states - two heavy and stable, 
opposite charge particles produced centrally plus 
two forward scattered protons. Using the phe- 
nomenological spectrum of [13], one finds the 
cross section for the two-photon production 
of about 0.43 fb. 

A usual method in pp studies to detect such 
heavy lepton like pairs relies on the use of 
the dE/dx variable, or time of flight measure- 
ments, and results in a poor reconstruction us- 
ing calorimeters and muon chambers [15] . On 
the contrary, in the two-photon analysis, almost 
all the kinematic information is available, again 
thanks to the detection of the forward scattered 
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Figure 8. Photon-photon invariant mass distribu- 
tion after the integrated luminosity L = 100 fb _1 for 
pp(77 — > T^T^)pp with Pt(?i ) > 10 GeV require- 
ment. 



analyzing the stau angular distributions in their 
center-of-mass reference system. 

6. Conclusion 

Photon-photon interactions provides an addi- 
tional and complementary handle on the BSM 
theories in general, and on the supersymmetric 
low mass models in particular. As it has been 
shown, the main advantage of the two-photon 
production analyses lies in the extra kinematical 
information - as the total center of mass energy, 
the missing energy, or the pair boost - obtained 
using the very forward detectors (at 220 m and 
420 m) at the LHC, as well as in the cleanliness 
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of the event signature. 

Discovery of the left and right-handed selec- 
trons and smuons after L = 25 fb _1 for the light 
MSSM scenarios is expected and mass measure- 
ment of these sparticles as well as of the lightest 
SUSY particle with a few GeV resolution should 
be possible. This assumes high resolution of the 
VFDs, also during high luminosity runs. 
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